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ABSTRACT 


As the costs for maintaining computer communication networks are rapidly rising, 
it is particularly important to design the network efficiently. The objective of this 
thesis is to model the minimum cost design of digital communication networks and 
propose a heuristical solution approach to the formulated model. 

The minimum cost design has been modeled as a zero-one integer programming 
problem. The Lagrangian relaxation method and subgradient optimization procedure 
have been used to find reasonably good feasible solutions. 

Although the reliability requirement for computer communication networks is as 
important as the cost factor, only the cost factor is considered in the context of this 


thesis. 
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I. INTRODUCTION 


Computer communication networks are becomming a vital part of today’s indus- 
trial, govermental, financial, and service organizations. The efficiency and quality of 
their internal operations and external services depends upon how effective their com- 
puter communication networks are. Design of computer communication networks 
requires the consideration of a number of factors such as performance constraints, 
reliability constraints, and costs. As the costs for maintaining computer communi- 
cation networks are rapidly rising, it is particularly important to design the network 
efficiently: i.e., the emphasis on keeping the costs Jow in the design of networks is 
becoming evident. 

Due to the inherent complexity of network design problems, the comprehensive 
optimal design of efficient networks is usually infeasible except for very simple net- 
works. The objective of this thesis is to model the minimum cost design of digiial 
communication networks and propose a heuristical solution approach to the formu- 
lated mocel. 

It is assumed that node locations are fixed by the required placement of terminals 
and network facilities. Therefore network design is predominately a question of where 
to put links and what their data carrying capacities should be. For a given set of n 
nodes, there are n(n-1)/2 possible links. It is seldom necessary for them to be fully 
connected. Consequently, there is an inherent design tradeoff between equipment 
simplicity and low cost with fewer links. Although the reliability requirement for 
computer communication networks is as important as the cost factor, only the cost 
factor is considered in the context of this thesis. 

The minimum cost design has been modeled as a zero-one integer programming 
problem. The Lagrangian relaxation method and subgradient optimization procedure 
have been used to find reasonably good feasible solutions. 

The organization of this thesis is as follows: Chapter II is an overview of digital 


transmission facilities and describes T1, T3 circuits and digital cross-connect systems. 





Chapter III formulates a model of the minimum cost design problem and introduces 
a Lagrangian relaxation method along with a subgradient optimization procedure. 
How to implement this technique in network design is also described. Chapter IV 
presents the results of computational tests. Finally, conclusion and recommendations 


for further research are offered in Chapter V. 





Il. DIGITAL TRANSMISSION FACILITIES 


A. Overview of Digital Transmission 

Digital transmission techniques were introduced in the Bell Svstem in the early 
1960s for efficient transport of voice signals. At that time, data for computer com- 
munications accounted for a small percentage of network traffic. Therefore, the basic 
digital transmission structure was centered around human voice communication. 

Although a plot of amplitude versus time for a speech sample shows significant 
high-frequency components, a successful telephone conversation requires less than 
4KHz of audio bandwidth. After being processed by a low-pass filter, a speech sampie 
loses its high-frequency components and is ready to be digitized. 

Since the maximum analog frequency to be reproduced is 4K Hz, an 8Kbit/s sam- 
pling rate was adopted according to the Nvquist theorem. By sampling the voice 
signal every 125 microseconds, its essential information is extracted from the ana- 
log format and is ready for digital encoding. This pulse-amplitude modulated(PAI) 
signal contains all the information in the original signal up to approximately 4K Hz. 

The modulation scheme that was chosen for the early digital transmission stan- 
dards is the easiest to implement. The pulse code modulation(PCM) is currently 
used worldwide for a standard voice communication. Figure 1 shows how the voice 
signal is encoded to digital signal. The pulse amplitude modulation(PAM) signal is 
quantized by mapping into discrete amplitude levels, each with a unique binarv code. 
Naturally, additional discrete mappings reduce the quantization error and improve 
the fidelity of the coded signal. The analog voice signal is sampled 8000 times per 
second with each sample measuring the amplitude of the signal, and sampled ampli- 
tude is encoded into an &-bit digital code. By creating 8&-bit codes 8000 times per 
second, the rate of a digitized voice signal is § multiplied by 8000, or 64Kbit/s. This 
forms the basic 64Kbit/s channel that is the foundation of the digital network. Such 


a digital channel is often referred to as a DS-O(digital signal level 0). 
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Figure 1: Pulse Code Modulation 





The original deploymant of digital transmission was driven by a desire to conserve 
copper pairs outside telephone offices. Each analog voice signal in the existing tele- 
phone network consumed a physical pair of copper wires from a subscriber location 
to a telephone central office. Once the basic block was digital, it became feasible to 
multiplex the digital signals together into a higher-order digital signal. The time- 
division multiplexing led to defining the next level of the digital hierarchy as being 
equivalent to 24 DS-0 signals. This level is referred to as DS-1 or T1. [Ref.13] 

The multiplexing function that created a DS-1 signal was originally handled by 3 
network element known as a digital channel bank or D-bank. Numerous devices, of 


course, offer T1 interfaces at these days. 


B. T1 System 
1. What is T1? 
The T1 system has become one of the most widely used high-capacity sys- 


tems for transmission of voice and data. Originally conceived as voice transmission 





technology in the early 1960s, it has evolved into a cost-effective and highly flexible 


means of transmitting both voice and data. 

Tl is based on time-division multiplexing (TDM) 24 users onto one phys- 
ical circuit. It usually involves a digital service unit(DSU) and a channel service 
unit(CSU) to accomplish the analog-to-digital conversion and the multiplexing. At 
the inception of T1, it was known that solid copper cable was capable of provid- 
ing frequencies about 100KHz. Consequently, the earlier frequency division multi- 
plexed(FDM) telephone system multiplexed 24 users onto one copper wire. Each 
user was assigned a different 4-KHz frequency band. The composite of 24-voice chan- 
nels totalled 96KHz, which was within the capacity of twisted pair wires. [Ref.11] 

The T1 system is designed around a 1.544Mbit/s rate, which was about the 
highest rate that could be supported across twisted pair for a distance of approxi- 
mately one mile. 

T1 has several advantages. The economic incentive to use T]1 can be enormous 
despite the apparent high initial cost cf equipment and the continuing costs for leased 
lines. An investment in a new network based on T1, however, can break even within 
one year. The end user organization operating a private Tl network enjoys the 
diagnostic power to identify and isolate faults and the abilitv to set up and change 
connections between points on the network at will. Another advantage is that Tl 
technology simplifies, and thereby improves corporate communications. Combining 
multiple voice and data channels over a single high-speed digital circuit eliminates 
many separate networks and many forms of network management. Therefore, T1 
makes network control more practical and improves reliability. With the availability 
of large amounts of bandwidth, under immediate control, the network operator is able 
to offer the corporation new services that were not practical in the past. Additional 


circuits are readily available through reconfiguration. [Ref.11] 


2. T1 Frame Format 
The majority of T1 offerings digitize the voice signal through PCM or adaptive 
differential pulse code modulation(ADPCM). Whatever the encoding technique used, 
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Figure 2: The Tl Based Frame 


once the analog images are translated to digital bit streams, many T1 systems then 
time-division multiplex voice and data together in 24 or 44 user slots within each 
frame. 

This section introduces the concept that the carriers transmit 24 voice or data 
channels together with TDM frames. The T1 carrier system provides multiplexing 
by sampling the 24 channels at a combined rate of 192.000 times per second. Figure 
2 shows how the 24 channels are multiplexed into a frame. The frame contains one 
sample from each channel, plus an additional bit for frame synchronization. Thus, 
the complete frame contains 193 bits (§ bits per channe] x 24 channels + ] sync. 
bit = 193bits). In addition. each of the 24 channels requires a 64Kbit/s rate (8000 
samples/sx § bit/sample = 64Kbit/s). 

The basic T1 frame consist of 24 eight-bit slots and a framing bit. The details 


of Tl frame format are characterized by the following concepts: 


Alternate Mark Inversion (AMI) 


ls density requirements 


Bit robbing 


e Binary 7 zero substitution (B7ZS) 
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Figure 3: Format errors Due to Bipolar Code Violation 
e Framing bit 


Alternate Mark Inversion 

The AMI concept is shown in Figure 3(a). The user may recognize this code as 
bipolar coding. The T1 system requires that a 1 pulse must be sent as an opposite 
polarity from the preceding 1 pulse, regardless of the number of Os in between the 
two ls. This bipolar code performs well because it has no direct current component 
and therefore can be coupled with transformers. The scheme also makes efficient use 
of bandwidth. 


1s Density Requirements 





Tl provides no separate clocking signal. The timing or clocking information is 
embedded in the data stream. At the receiving end, the clock is recovered from the 
data stream by the detection of 1 pulse. If the T1 data stream has insufficient 1 
pulses embedded, the receiver can no longer produce reliable timing out. 

To overcome this problem, a certain number of 1s must be present to ensure proper 
timing. This concept is called 1s density. The T1 facilities require that no more than 
15 consecutive Os shall be present in the frame. An additional requirement is that 
there must be at least three 1s in every 24 bits. 

Due to noise, mechanical failures, etc., the bits may become distorted, which can 
cause a violation of the AMI rule. Bipolar violations or excessive errors are known as 
format errors because the errors are not in conformance with the required T1 format. 
A bit distortion may not cause a format error. As Figure 3(b) through 3(e) show, the 
nature of the error determines if the bipolar violation is detected. The bipolar signal 
is altered with errors of omission or errors of commission. 

Bit Robbing 

The least significant bit of every time slot in the ’signalling frames’ is overwritten 
to encode the control signals. This concept is called bit robbing. 

For the transmission of data, the 8th bit is unavailable. Consequently, the majority 
of Tl and related systems use a 56Kbit/s transmission rate instead of the 64Kbit/s 
rate actually available. 

Binary 7 Zero Substitution 

To establish 1s density in the T1 signal, a technique called B7 zero code suppres- 
sion is used. As previously discussed, a 1 is substituted in a T1 frame to prevent an 
occurrence of more than 15 consecutive Qs. It is possible that this substitution can 
occur with any bit in the frame. A voice channel can tolerate the loss of the 7th bit 
to B7 substitution. However, a data channel cannot operate with this arrangement. 

It is not necessary to rob bits in a data channel since there is no telephone signal. 
However, it is quite possible to have all 0s in a data channel. Therefore, the data 


could be corrupted by the B7 zero code suppression technique. The insertion does 
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Figure 4: T1 Facility Termination 


not affect the voice channel, but it certainly affects a data channel. To prevent this 
problem, a data channel generally only uses 7 bits and any 1 substitution can be 
placed into the 8th bit. This prevents the B7 zero code suppression from corrupting 
the data channel. 
Framing Bit 

To decode the incoming data stream, a receiver must be able to associate each 
sample with the proper channel. At a minimum, the begining and ending of the frame 
must be recognized. The function of the framing bit provides this requirement. This 
bit is located in the 193rd bit of each frame. It is not part of user’s information, 
but added by the system for framing. The use of framing bit varies significantly 


depending on the type of Tl and of the technology. [Ref.17] 


3. T1 Networks 
(A) T1 Interface 
Figure 4 shows T1 facility termination. The T1 facility termination enters 
the users premises as two wire pairs for data. There may be additional wires for 


loopback control and testing of unframed lines. The channel service unit (CSU) 
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Figure 5: Digital Cross-Connect System 


provides loopbacks to the network and is the last regeneration point for incoming 
signals. DSX-1 as an interface is well defined in terms of electrical and mechanical 
features. A digital service unit (DSU) matches various terminal interfaces to the 
DSX-1 standard and usually is part of the user’s terminal equipment: multiplexer, 


bandwidth manager, PBX, computer, etc.[Ref.11] 


(B) Digital Cross-Connect Systems 
Digital cross-connect systems (DCS) are electronic switches able to take 
bit streams out of the Tl-carrier and route them to a given output line. DCS allow 
remote access to and cross-connection of large bandwidth circuits. Figure 5 shows 
DCS. DCS can be divided into three primary subsystems. 
At the core of the DCS is the Matrix Subsystem where the actual chan- 


nel level (DS-0) cross-connections take place. DCS systems are controlled by the 
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Figure 6: Data Bypass 


Administrative Subsystem, where high-level processors manage the overall system 
functions. Other common elements such as system synchronization and user inter- 
face/control port access reside within this subsystem. Finally, there is the Interface 
Subsystem, which terminates DS-1 or higher facilities. Interface subsystems are orga- 
nized into growth increments, or units, which accommodate anywhere from 28 to 32 
slot positions each. Aside from facility terminations, these slot positions can be used 
to accommodate features such as sub-rate date multiplexing and bridging of digital 
services. 

Figure 6 shows data bypass based on DCS. Data bypass resembles the 
drop and insert function of the voice channel bank. What distinguishes the two is 
that data bypass offers variable bandwidth allocation and additional flexibility in 


routing connection among many more T1 links.[Ref.11] 


(C) T1 Networking 
Lower cost for digital lines and for customer premises equipment have 
made all-digital networking not only possible, but attractive. Every location is on 


the digital backbone, through smal] sites may be served by 56Kbit/s or other data 
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Figure 7: T1 Networking 


link speed less than a full T1. 

The T1 networking are shown in Figure 7. There are two main benefits 
in Tl networking. Trends in tariffed costs are up for analog facilities and down for 
digital facilities. As time goes on, it takes fewer analog lines to justify a digital line. 
Currently, as few as 3 or 4 leased voice grade lines cost as much as a 56kbit/s digital 
service. As few as two or three 56K lines cost as much as a T1]. The result is that we 
can cost-justify digital service in more and more locations. Although there exist var- 
ious types of nodes with a common network management system, everything can be 
brought under one network management system because T] networking is all-digital. 


[Ref.11] 


i2 


C. T3 SYSTEM 
1. What is T3? 

Many corporate managers are administering T1 networks that are growing 
rapidly with new users asking for more and more bandwidth. Increasingly, data 
communications managers are eyeing T3 network as a solution to the problems raised 
by this expansion. 

T3 refers to a digital transmission facility running at DS-3 data rate of 
44.736Mbit/s. It is equivalent of 28 Tls, and is more economical than a much less 
number of Tls. For a distance of less than 50 miles, only four T1ls are required to 
break even. 

T3 circuits usually use microwave and optical fiber as transmission media. 
Optical fiber has a much higher transmission capacity than digital microwave. Dig- 
ital microwave has some problems. First, the electromagnetic spectrum is strictly 
licensed by the Federal Communications Commission. Second, strict spectrum al- 
locations imply that digital microwave cannot be readily upgraded in bandwidth. 
While a 45Mbit/s optical fiber can be readily converted to 565Mbit/s, a 45Mbit/s 
digital microwave setup is probably destined to remain at 45Mbit/s. Finally, digital 
microwave devices are inherently limited to line-of-sight distance. 

Many T3 networks mix optical fiber and digital microwave. The optical fiber 
can be used in high-density areas, with microwave hops to lower-density remote lo- 
cations. Alternatively, the main network can be placed entirely on fiber, with T3 
microwave used for emergency restoration in the event of a cable cut. Dual-media 
transmission is especially attractive to service-critica]l networks such as banks, airlines, 
and other on-line transaction-processing companies. 

T3 was originally designed to be used as interexchange trunk within the public 


telephone network.[Ref.13} 


2. T3 Format 
By the end of the 1970s, a way around the limitations of the copper pair 


was on the horizon. Transmission technologies based on the nearly limitless band- 
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width of optical fiber were emerging from the laboratory, and a new layer of the 


digital hierarchy seemed appropriate. Based on technology available at the time, 
an asynchronous DS-3 rate was defined as the combination of seven DS-2 signals of 
6.176Mbit/s. Framing and stuffing bits, as well as rudimentary error checking and 
internal communications, were added to the DS-2 tributaries to generate an aggregate 
bit rate of 44.736Mbit/s. 

Notice that the information content, or payload, of a DS-3 is equal to 672 
times 64Kbit/s = 43.008Mbit/s. The additional 1.728Mbit/s of overhead represents 
the sum of the DS-1 framing bit, the DS-2 framing and stuffing bits, and the variety of 
DS-3 overhead bits. In terms of network efficieny, the DS-3 format devotes 96 percent 
of the transmission bandwidth to payload, with approximately 4 percent overhead. 
[Re f.13] 

Although the DS-3 was designed to be created from multiple DS-2 signals, 
there was an obvious inefficiency in using two separate devices for the DS-1-to-DS-2 
and DS-2-to-DS-3 multiplexing functions. 

There are few choices when it comes to T3 formats. The long established 
standard in the public network is now called M13. But for technical reasons there 


are new proposals for formats called Syntran and SONET. 


(A) M13 Format 

The name M13 comes from the fact that this multiplexer connects the 
DS-1 to the DS-3 level. One DS-3 contains 28 DS-1s. They are combined in two steps. 
The M13 process first joins four DS-1 lines into DS-2 bit stream at 6.312Mbit/s. This 
step is more than a straight time division multiplexer function. M12 bit-interleaves 
the four inputs, but it also adds 136Kbit/s of overhead and justification or bit stuffing. 

The second step within M13 combines seven DS-2 streams into one DS-3. 
Like M12, the M23 process is bit-interleaved, with justification. That is, the DS-3 
clock need not be synchronized with any DS-2 or DS-1 line speed. M13 overruns the 
DS-2 clock, and bit stuffs or adds extra bits to fill out the DS-3 rate of 44.736Mbit/s. 


Whatever the original reasons for the two-step multiplexing, it creats a 
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problem for us today. There is no reliable way to tell where a given DS-0 is located 
within the DS-3 bit stream. This makes it very difficult to extract just one DS-0 
because it will have no definite location within the DS-3 frame. 

To get at an individual channel you must typically reverse the two step 
and demux the DS-3 down to individual DS-is. The need for a DS-3 digital access 
and cross-connect system(DACS) is clear. To get simple switching capability among 


DS-3 lines, a better format is required. [Ref.11] 


(B) Syntran 
The Syntran has been proposed by Bellcore. Syntran describes how to 
multiplex DS-0 and DS-1 signals into a DS-3, in a fully synchronous format. 
The resulting bit stream has the same rate as the standard. Syntran offers 


two important features. 


e Synchronous format allows easy DACS switching of DS-0 and DS-1 channels 
among DS-3 ports. 


e The aggregate bit rate is compatible with the large number of existing DS-3 


transmission facilities installed in North America. 


Unfortunately for the global network, DS-3 is seldom found outside of 
North America. And while M13 and Syntran are defined as electrical] interfaces, the 


new installations of facilities near that speed are mostly optical fiber.[Ref.13} 


(C) SONET 

SONET(Synchronous Optical Network) is much more than DS-3. The 
goal is a standard for interconnection of different national networks at speeds from 
about 150Mbit/s up to several gigabits/s. The DS-3 level is only one portion of the 
standard, intended for network access from customer premises. 

SONET operates at multiples of T3 bandwidth. Initially, products will be 
offered at the following bit rates (OC stands for Optical Carrier): OC-1, 51.84Mbit/s; 
OC-3, 155.52Mbit/s; OC-12, 622.08Mbit/s; OC-48, 2.49Gbit/s. 
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As asynchronous standard, SONET is well-suited for switching tributary 
signals within a higher bandwidth pipe. Initially, switching will be limited to T1 and 
DS-0 signals, but other service offerings will be defined in the future. 

SONET is fully backward-compatible with the ANSI synchronous T3 for- 
mat; it is partially backward-compatible with order asynchronous T3 equipment such 
as M13s. It does not affect users of digital microwave radio, since there is little need 


for a synchronous standard higher than T3 for these applications. [Re f.13] 


3. T3 Networks 

Implementing a T3 corporate backbone provides users with a measure of flex- 
ibility and control over the network. At the low-speed end, circuits may be either 
56Kbit/s or 64Kbit/s[Ref.13]. These may be used for voice or data services. 

The capability to switch the variety of circuits allows the user to put up and 
take down circuits of differing rates on the T3 backbone as requirements change, 
without having to work through an external network provider. This can take the 
form of reducing lead time for a circuit order, time-of-day circuit changes, or quick 
response to a temporary overload condition in some part of the network. 


T3 networks are based on one or more of the following types of DCS’s. 


1. Point to point terminals: M13 multiplexers with 28 T1s in and one DS-3 


out. 


2. Multiple M13 multiplexers in a single chassis: something like a 48-port D4 
channel bank that followed the 24-port D3. 


3. Networking of multiple DS-3 lines into one hub with the ability to cross- 
connect DS-1s and even DS-0s within that hub. 


Level 3 generally is possible only with one of the synchronous formats, not 
with classic M13. Devices with this performance are installed at operating companies. 
Network topologies for T3 backbones will evolve from point to point up to 


the partial mesh or interlocking rings familiar from T1 nets. There will be many 


16 





fewer sites that can cost-justify multiple DS-3s than there wil: be Tl hubs. The 


combination likely to be most common is a sparse DS-3 backbone surrounded by a 


T1 network. [Ref.11] 





III. MODEL FORMULATION AND SOLUTION 


This chapter discusses how to formulate the minimum cost design of T1/T3 net- 
works, and explains the heuristic procedure for optimal solutions based on the La- 


grangian relaxation method. 


A. MATHEMATICAL FORMULATION 

In order to model the minimum cost design of a network, we consider several 
things. The amount of traffic between each node and every other node is given and 
we define the set of circuit types each of which is a mix of T1, T3 and fractional T1. Ifa 
network ccntains n nodes, there would be n(n-2)/2 communicating source-destination 
pairs that select routes from sets containing from a few thousand to hundreds of 
thousands of candidate routes, leading to very large problem sizes. Therefore, to 
simplify the computation, we assume that the maximum length of the path for each 
pair is less than a specified number of hops. For example, the network designer 
may specify that routes between two nodes should pass through no more than two 
additional nodes. 

We introduce the following notations to present the mathematical formulation of 


the problem. 
P : The index set of the communicating source-destination pairs 


R, : The index set of candidate routes for the source-destination pair p € P. The 
set of candidate routes is generated by a route generation algorithm. A route 


is characterized by the ordered set of links from source to destination node. 
L : The index set of candidate links in the T1/T3 network. 
T : The index set of circuit types. 


X pr : A decision variable which is one if the source-destination pair p takes the route 


r and zero otherwise. 
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Zu : A decision variable which is one if the link / is supported by the circuit type t 


and zero otherwise. 
Q, : The number of DS-0 channels that can be carried over a circuit type t. 
a, : The number of DS-0 channels demanded by the source-destination pair p. 
Ci, : The costs for leasing a circuit type t on link J. 


I,, : An indicator function which is one if link / is used in route r and zero otherwise. 


Using the above notations, we formulate the minimum cost design problem as 


lollows: 
Problem P: 
Zp =min)_ >> CuZu (1) 
leL teT 
subject to 
De Agee a Vp Ee P (2) 
r€Rp 
Ss oD GolriXge < we Vl = L (3) 
p€P reR, teT 
YZ <1 Vie L (4) 
teT 
0< Xy <1, 0< Zn <1 integer (5) 


The constraint (2) ensures that each source-destination pair selects exactiy one 
route for transmission. The constraint (3) ensures that the number of DS-0 channels 
on each link does not exceed its capacity. The selection of at most one circuit type 
for each link is ensured by constraint (4). Problem P is a linear zero-one integer 
programming problem which can be shown to be NP-hard. 

Based on the current status of the computational ability, the idea of solving 
Problem P optimally is rejected. The idea of using piecewise linearization of the 


objective function and solving the linear programming problem is also rejected for 
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the reason that the linear relaxation of Problem P is very large and its optimal solu- 


tion requires a significant computational effort and may generate a fractional solution. 


B. LAGRANGIAN RELAXATION 

_we of the most computationally useful ideas of the 1970s is the observation that 
many hard problems can be viewed as easy problems complicated by a relatively small 
set of side constraints. Dualizing the side constraints produces a Lagrangian problem 
that is easy to solve and whose optimal value is a lower bound for minimization 
problems on the optimal value of the original problem. The Lagrangian problem 
can thus be used in place of a linear programming relaxation to provide bounds in a 
branch and bound algorithm. The Lagrangian approach offers a number of important 
advantages over linear relaxation. 

Lagrangian methods had gained considerable currency by the 1974 when Geoffrion 
coined the perfect name for this approach - Lagrangian relaration. Since then the 
list of applications of Lagrangian relaxation has grown to include over a dozen of the 
most infamous combinatorial optimization problems. For most of these problems, 
Lagrangian relaxation has provided the best existing algorithm and has enabled the 
solution of problems of practical size. [Ref.3] 

In the previous section, we already made a mathematical] formulation for minimum 
cost design problem. We know that there are two natural Lagrangian relaxation for 
Problem P. The first is obtained by multipling the constraint (3) by a vector of 


positive Lagrange multipliers, U;, | € L, and adding them to the objective function. 


Problem D,: 
Zni(U) = min > 52, CuZit + ie U; > > ap],1X pr = Ds O21 
leL t€T leL pePreR, t€T 
= min » (Cu —UiQ:)Zin + is > >. Uap] yi X pr (6) 
leL t€T lEeL peP reRp 


subject to 








he, Ae ae Vp e P 
reERp 
>> Zu ae | Vie L 
teT 


0< X, <1, 0< 2, <1 integer 


The set of feasible solutions for Problem P is a subset of the set of feasible solutions 
for the Lagrangian relaxation D,. The positivity of U ensures that in any feasible 


solution for Problem P, the expression 


a U; De > pl yj X pr — > O:Zi 
teT 


leL pE€Pré€Rp 


is nonpositive, and thus, the value of the objective function in (6) is never greater 
than the value of the objective function in Problem P. Thus, whenever Problem 
P has a feasible solution, Zp,(U) is less than or equal to Zp. For each vector of 
multipliers U, Zp,(U) is a lower bound for Zp. 

The second relaxation is obtained by dualizing constraint (2) with nonpositive 


Lagrange multipliers,V,, p € P. 
Problem D2: 
Zp(V) = min) >> CuZu + a Vy b> Age = | (7) 
l€L t€T peP reéRp 


subject to 


ys > Onli Xpr < wey 77, VIE L 
p€P r€Rp teT 
Zu <1 vie L 
teT 


0< Xp <1, 0< Zu <1 integer 


Comparing above two Lagrangian relaxations, we see that the second relaxation is 


harder to solve but might provide better bounds. This thesis uses the first Lagrangian 
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relaxation for implementation and computational results. 


C. SUBGRADIENT OPTIMIZATION PROCEDURE 

Several methods have been suggested in the literature for computing a vector of 
optimal Lagrange multipliers. These include various versions of the simplex method 
implemented using column generation procedures, dual ascent and multiplier ad- 
justment procedures, and a subgradient optimization procedure[Ref.3]. Subgradient 
optimization methods have been shown to be very effective in a variety of combinato- 
rial optimization problems such as the traveling salesman problem, topological design 
of computer communication systems, and lot-sizing algorithms for complex product 
structures|Ref.1]. Therefore, we decided to employ subgradient optimization tech- 
niques. 

The subgradient method is a brazen adaptation of the gradient method in which 
gradients are replaced by subgradients. 

Let (Xp-(U), Zi(U)) be the optimal solution to the Lagrangian problem D, for a 
fixed vector U. A subgradient is the vector with coordinates 


n(U) = > YS apliXp(U)— 32 Q:Z(U) WEL. 


p€P reéR, teT 


Given an initial value U°, a sequence U* is generated by the iterative formula 
k k k 
U;t" = UF + ten 


where t, is a positive scalar stepsize. It is known that Us converges to the optimal 
Lagrange multiplier” provided that t, converges to 0 and > ¢; diverges. The stepsize 
most commonly used in practice is 
a Xk [Zp 7 Zn:(U*)| 

Ml (I? 
where i, is a scalar satisfying 0 < A, < 2 and Zp is an upper bound on Zpi(U), 
frequently obtained by applying a heuristic to Problem P. Often the sequence , is 


determined by setting Ao = 2 and halving 4, whenever Zp,(U) has failed to increase 
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in some fixed number of iterations. This rule has performed well empirically, even 
though it is not guaranteed to satisfy the sufficient condition given above for optimal 
convergence. 


The steps involved in the subgradient optimization procedure are as follows. 
1. Initialization: 
(a) Using a heuristic, compute an overestimate Zp of Zp or set Zp to an 
arbitrarily large value. 


(b) Select an initial set of multipliers U°; and set K, the iteration counter, to 
0, and the improvement counter to 0; and \ a parameter for adjusting the 


stepsize to \° (an arbitrary positive initial value, e.g., 2). 
2. Solving the Lagrangian problem: 


(a) Increment the improvement and iteration counters by 1. 


(b) Solve the Lagrangian problem using U* as the Lagrange multipliers. Thus, 
obtain Zp,;(U*) and X*, Zi. 


prs 


3. Testing and updating parameters: 


(a) If Zp,(U*) is greater than the current best value of Zp,(U) then replace 
the current best value of Zp,(U) by Zp, (U*) and set U* to U*. Also reset 


the improvement counter to -1. 


(b) If X‘ and Zj, are feasible for Problem P, compute the value of its as- 
sociated objectivefunction for Problem P. If this value is less than the 


current value of Zp. then set Zp to this value. 


(c) If the improvement counter has reached a prespecified upper limit, then 
set \ to \/2, &* tow", the improvement counter to 0 and go to 2(a). 

(d) If the iteration counter has exceeded a prespecified limit, if 4 is less than 
a prespecified limit, if t, is less than a prespecified limit, or if (Zp - 


Zp1(U*))/Zp,(U") is less than a prespecified error tolerance, then stop. 


23 


4. Updating the multipliers : 


(a) Compute a new subgradient. 
of => YS aplnXt-SQZi Viel 
p€P reRp t€T 
(b) Compute the new stepsiz. 
i [Zp - Zp. (U*)] 
; oF IP 


(c) Compute the new multipliers. 
Ur =UFtigt =e 
(d) set k tok +1. 


5. Go to 2. 


D. SOLVING THE LAGRANGIAN PROBLEM 

The computational efficiency of the subgradient optimization procedure depends 
on our ability to solve efficiently the Lagrangian problem which is generated in step 
2(b) of the subgradient optimization procedure. Fortunately, for a fixed set of multi- 
pliers, the Lagrangian problem is separable into subproblems which are readily solved. 

Since there are no coupling constraints between the X,,, and Z, variables, the 
Lagrangian problem D, can be decomposed into two subproblems as Zp,(U) = 
Zpu(U) + Zpi2(U) where 


Subproblem Dy,: 


Zou(U) = min >> > So Ural X pr (8) 
l€L p€P r€Rp 
subject to 
> Ap = 1 VpeP 
rE€Rp 


0< Xp, <1 integer 
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and 


Subproblem Dy2: 


Zpi2(U) = min > Yo (Cr — UiQs) Zit (9) 
leL teT 
subject to 
Soe oO WEL 
teT 


O< Zn <1 integer 


Subproblem D), can be separated into | P | subproblems, one for each pair, where 


the subproblem associated with the p* pair is 


Z5u(U) = mina, >> pu Xo 


rER, LIEL 


subject to 


> Xr = 1 
réeRp 
0< Xp, < linteger 
and then Zpu(U) = opeP 2513 (U). 
Subproblem D2 can also be separated into | L | subproblems, one for each link, 


where the subproblem associated with the /"* link is 


Zpi2(U) = min (Cu — UiQ1) Zn 
teT 


subject to 


S- Zu <tl 


t€T 
0< 2Zn<1 integer 
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and then Zp12(U) = Der Z2(U). 


E. IMPLEMENTATION PROCEDURE 

To test the applicability of the subgradient optimization procedure to network 
design problems and to obtain estimates on the quality of bounds and solutions 
generated, this thesis implemented the procedure proposed in previous section in 


turbo Pascal on IBM PC computer. 


1. System Flow 

The program consists of two main parts. The first is data generation, and the 
second is problem solving. Figure 8 shows the system flow. The program is initiated 
by entering the number of nodes. It then creates data file with source-destination 
node pairs, and their distance and amount of data traffic. The user is then prompted 
to enter a number that he/she wantsa to generate a set of candidate routes. Next, 
the user is prompted to specify a number of subgradient parameters. 

Two types of iterations, major and minor, are defined in the program. A 
major iteration consists of many minor iterations, each of which is an evaluation 
of the current value of the objective function, a choice of subgradient direction and 
modification of the value of U. The minor iterations terminate when one of the 
stopping conditions is satisfied such as the difference between the best feasible solution 
and the lower bound based on the Lagrangian dual is small, the subgradient is smal] 
or an upper limit on the number of iterations has been reached. 

Throughout the optimization process, the program updates the Lagrangian 
value, the multipliers generated in each major iterations and the best feasible solution 
yet generated. 

Based on the output generated after a major iteration, the user may terminate 
or reenter the subgradient optimization procedure by changing the number of routes, 
iteration counter or improvement counter. This iterative process continues until the 


solution does not get improved. 
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MAIN BODY 


GENERATE DATA SOLVE PROBLEM REPORT 


INPUT DATA READ DATA 


MAKE PAIRS INI. SOL. 
MAKE CIRCUIT INI. SUB PARA 
GET SOLUTION 





Figure 8: System Chart 











2. Heuristic 

The subgradient optimization procedure requires an overestimate Zip on 
Zp1(U*) used to update the Lagrange multipliers. This overestimate can be a feasible 
solution generated by a heuristic. Figure 9 shows the heuristic that has been used in 
this thesis to generate good feasible solutions. 

In each minor iteration of the subgradient optimization procedure, dual solu- 
tion X* and Z;; are obtained. These solutions might be infeasible because a set of 
capacity constraint (3) was relaxed. This constraint, therefore, may be violated by 


the dual solution Xp and vA 
Yo DY ples X pe > Do Ze Je L. 
peP réeRp téT 
By correcting this infeasibility situation for all 1 which has been violated, 
2 a Ap] 1X pr < yO Zit VIEL, 
peP reR, t€T 


the solution becomes feasible. 
If this feasible solution is smaller than previous best feasible solution which 
is upper bound, it replaces the previous solution and becomes the new best feasible 


solution. 
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INPUT : Infeasible Dual Solution X*. and Zi 
OUTPUT: Feasible Solution 


BEGIN 
OBJECTIVE := Q; 
FOR each Link ! 
BEGIN 
t := 0; 
WHILE ap],1Xpr > Q: DO 
t:=t+1; 
Zi := t; 
OBJECTIVE := OBJECTIVE + CuZu 
END; ( FOR LOOP ) 
IF OBJECTIVE < UPPER-BOUND THEN 
UPPER-BOUND := OBJECTIVE; 
END. ( HEURISTIC ) 


Figure 9: Heuristic 








IV. COMPUTATIONAL RESULTS 


The solutiop approach proposed in Chapter III has been tested using network 
design problems. 


A. Computational Experiments 

The procedure was tested under two situations: Under one situation [case A], the 
cost of T3 circuit was 9 times of cost of T1. The other situation [case B] was such that 
the cost of T3 circuit was 5 times of Tl. Under each situation, we have tested several 
networks that have different number of nodes. To find better feasible solutions, we 
have tested this procedure several times changing the number of candidate routes — 
8, 12, 16 and 20 candidate routes. 

The distance and amount of traffic for each of source-destination node pairs and 
circuits information were randomly generated by the random number generation pro- 
cedure. Figure 10 shows one of sample data used. 

The procedure has the number of subgradient parameters. For the test, we ini- 
tially set \° to 0.5 and U® to 0. The stopping criteria parameters are set as follows: 
stepsize limit to 0.0001, error limit to 0.005, and iteration limit to 100. 


B. Summary of Computational Results 

The results of the computational tests on the five randomly generated problems 
are summarized in the Table 1 and 2. The values shown in the tables represent the 
feasible solutions for network with 6, 10, 15, 20 and 25 nodes. The asterisk means 
that the previous feasible solution is not improved. 

While testing the procedure, we had some problems with the size of main memory 
and the processing times. One of the problems was memory overflow error when we 
were trying to enlarge the number of node, the number of candidate route of each pair, . 


and/or the number of circuit. Since it takes O(n"~*) times to locate all candidate 
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node pair information 
source dest. distance traffic 


1 2 443 0 
1 3 592 60 
1 4 488 71 
1 5 516 0 
1 6 648 4 
2 3 149 .¢) 
2 4 607 0 
2 5 319 53 
2 6 539 0 
3 4 695 0 
3 5 351 61 
3 6 562 70 
4 5 371 64 
4 6 298 46 
5 6 220 97 


circuits info. 
capa. cost 


6 10 

12 25 
18 19 
24 21 
30 31 
1554 625 
1560 627 
1566 637 
1572 642 
1578 646 


Figure 10: Sample Input Data 
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Time (sec) 


oes for 20 candidate routes 





for 10 candidate routes 








12 14 16 
No. of nodes 


10 18 20 ----:- 


Figure 11: Graph of Time Complexity 


routes, we limited the maximum number of hops of candidate routes to 4 and the 
maximum number of nodes to 25. Figure 11 shows how much time is required as the 
number of nodes increases. 

Table 1 and 2 summarized the results obtained for networks with the different 
candidate routes. In most of cases, the feasible solution improve when the number of 


candidate routes increases. 
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N, : Number of nodes 

N, : Number of candidate routes 
for each node pair 

* : Indication for nonimprovement 


Table 1: Summary of Feasible Solutions : Case A 
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Table 2: Summary of Feasible Solutions : Case B 
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V. CONCLUSION 


We have made a mathematical formulation for minimum cost design of digital com- 
munication network. This thesis used the Lagrangian relaxation and the subgradient 
optimization procedure to obtain lower bounds and feasible solutions in minimizing 
the cost of network design with given node location and traffic of each node pairs. 
Computational experience indicates that this procedure is effective in most instances. 
The quality of bound generated depends to a large extent on the limits set on the 
improvement and iteration counters. Low settings terminate the procedure before 
it reaches a good feasible solution, while high settings consume excessive computing 
resources. Good settings of these parameters require careful balancing between these 
two consideration. 

The difference between upper bound and lower bound was fifty to seventy percent 
in our experiments. Our conjecture is that the other constraint relaxation is required 
to get a better bound. This thesis did not explore that constraint relaxation. 

Our model can be extended by considering reliability constraints, which is left for 


future research. 
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APPENDIX A: SOURCE CODE 


A. DATA STRUCTURE DESCRIPTION 


{ REKREKKEEKEKEKEKREKKEEKEEEEEEKREERERERREEREKEREKEKEREEEREKEEE } 


{* TITLE : Minimum Cost Design of Network *)} 
{* AUTHOR : Hwang, Yong Goo *} 
{* DATE : 13 June 1990 *} 
{* REVISE : 15 Nov. 1990 *) 
{* DESCRIPTION : This program is to get best *} 
(* feasible solution for minimum cost design of *) 
{* network. It has been implemented in Turbo *} 
{* Pascal and run on IBM PC computer. *} 
{* System flow is as follows: *) 
{* 1. Generate data and store it to ’net.dat’ *} 
{* file. *) 
{* - node pairs information *)} 
{* - circuit information *} 
{* 2. Get data from ’net.dat’ file *} 
{* - make route table *} 
{* 3. Solve the problem *)} 
{* 4. Out report report +) 


{ KHKKKKKKKKKKEKKKKEKKKEKKKKEKEKEKEKKREKEKKEREKEKEEKEKEKEKKEKEE | 


PROGRAM NETWORK_DESIGN; 


USES CRT, GRAPH; 


CONST 
VLARGE = 1€20; 
MAX_I = 20; { maximum No. of nodes } 
MAX _P = 190; 
MAX_R = 6000; 
MAX T = 200; { maximum size of circuit } 
max_ne = 20; { maximum No. of candidate routes } 
F_NAME = ‘net.dat’; 
FILE_NAME = ‘out.dat’; 
TYPE 
NODE = 1..MAX_I; 
PAIR = 0..MAX_P; 
CNFG = 0..MAX T; 


POINT = ARRAY[1..100] OF INTEGER; 
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CIRCUTE = RECORD 
FT_1 : INTEGER; 
T 1 : INTEGER; 
T 3 : INTEGER; 
CAPACITY: INTEGER; 
COST : INTEGER; 
END; 


CIRCUTE_TABLE = ARRAY[0O..MAX_T] OF CIRCUTE; 


SD_PAIR = RECORD 


S : INTEGER; 
D : INTEGER 
END; 


PAIR_TABLE = ARRAY[1..MAX_P] OF SD_PAIR; 


VAR 


CNFG; 


Be 
p: 
t: 
r,ni,np,nr,nt,nc : 


INTEGER; 


ORGN,DSTN: array[{1..MAX_P] of NODE; 
DISTANCE, TRAFFIC: array[1..MAX_P] of REAL; 
CAPA,RATE: array[0..MAX_T] of REAL; 


ROUTE: 
px, ax: 
pz,az: 


ITRCTR, ITRLIM, IMPCTR, IMPLIM: 


array(0..MAX_R,1..4] of PAIR; 
array[(1..MAX_P] of 0..MAX_R; 
array[(0..MAX_P] of CNFG; 


INTEGER; 


UPPBND, LOWBND, pobj,dobj: REAL; 
NORM, ERRORLIM, LAMBDA, LAMBDALIM, STEPSIZE: REAL; 


ALPHA, BEST_ALPHA, SUBGRAD: 
OPTIMAL, TERMINATED, FRAG1: 
DATA_FILE : TEXT; 
POINT X, POINT_Y : POINT; 
T TABLE : CIRCUTE_TABLE; 
S_D_PAIR : PAIR_TABLE; 
ANS1 : CHAR; 


array[0..MAX_P} of REAL; 
BOOLEAN ; 
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B. DATA GENERATION PROCEDURE 


{ RERRKKAEEKEKKEEEKEREKREEEREEKEREEKEREEEREEEEEKEKEKREKKEKE ) 


{* TITLE : Generate Data *)} 
{* AUTHOR : Hwang, Yong Goo *)} 
{* DATE : 15 June 1990 *} 
{* REVISE : 15 July 1990 *)} 
{* DESCRIPTION : This procedure will make a data *) 
{* file that consists of node pair informations *} 
{* and circuit informations with random genera- *} 
{* tion procedure *) 


{ RHKEKKRAKEKKEKEKEAEKKKKEEKEKEKEKEKEEKEEREEEEKEKEEKEKEKKKKE } 


PROCEDURE GENERATE_DATA; 


VAR 
t : INTEGER; 


{ RRARRRRERAERERERERRERRERERREKEEREREEEEEEKEKEKEE ) 
PROCEDURE T_TABLE_GENERATION; 


VAR 
i, j, k : INTEGER; 


k : 

EGIN 
T_TABLE(nt}].FT_1 := k; 
T_TABLE(nt}.T_1 := j; 
T_TABLE[nt).T_3 := i; 
T_TABLE{nt].CAPACITY := k*6 + j*24 + i*672; 
CASE k OF 


1 : T_TABLE[nt].COST := 10 + j*21 + i*219; 
2 : T_TABLE(nt}.COST := 15 + j*21 + i*219; 
3 : T_TABLE[nt].COST := 19 + 3*21 + i*219; 
ELSE 
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T_TABLE(nt].COST := j*21 + i*219; 
END; 
nt := nt + 1; 
END; 


END; 


{ REKRKKKAERREKERREREEEEEEREEEEKEREKKEKEREKEEEE ) 
FUNCTION RANDINT(MIN,MAX : INTEGER): INTEGER; 


{* The RANDINT function suppliers a random integer within 
the range of the specified arguments min and max. 
For example, RANDINT(10,20) gives a random integer 
between 10 and 20. *} 
CONST 
MININT = -5000; 


VAR 
RANDRANGE : INTEGER; 


BEGIN 


IF (MIN < MININT) THEN 

MIN := MININT; 
RANDRANGE := MAX - MIN + 1; 
RANDINT := RANDOM(RANDRANGE) + MIN; 


END; 
{ RERAEKRERERKKEEREREREREEEES ) 
PROCEDURE NODE_GENERATION; 


VAR 
n : INTEGER; 


BEGIN 


FOR n := 1 TO ni DO 
BEGIN 
POINT _X(n] := RANDINT(50,1000) ; 
POINT Y{n] := RANDINT(30,800); 
END; 
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END; 


{ RHERERREKKHEEKREKEKREKKERER } 
PROCEDURE PAIR_GENERATION; 


VAR 
3, k, m, p : INTEGER; 
DISTANCES : REAL; 
TRAFFIC_VOLUMN, TEMP_TRAFFIC_VOLUMN : INTEGER; 
X, ¥ : REAL; 


BEGIN 
p := 0; 
FOR j := 1 TO (ni - 1) DO 
BEGIN 
m:=j+ 13 
FOR k := m TO ni DO 
BEGIN 
Pt=p+ il; 
S_D_PAIR[p].S := j; 
S_D_PAIR[p]}.D := ki 
X := POINT XJ] - TOINT_X[(k]); 
Y := POINT _Y[(j] - POLIT _Y[k]; 
DISTANCES := SQRT((X * X) + (Y¥ * Y)); 
TEMP_TRAFFIC_VOLUMN := RANDINT(~30,100) ; 
IF (TEMP _| TRAFFIC |_VOLUMN < 0) THEN 
TRAFFIC_VOLUMN s= 0 
ELSE 
TRAFFIC_VOLUMN := TEMP_TRAFFIC_VOLUMN; 
WRITELN(DATA_FILE, 3:5, k:5, DISTANCES:10:0, 
TRAFFIC_VOLUMN:10) ; 
END; 
END; 
END; 


{* main part of data generation *)} 
BEGIN 


ASSIGN (DATA_FILE, F_NAME) ; 

REWRITE (DATA_FILE) ; 

RANDOMIZE; 

WRITELN (’ Enter number of NODEs? ’); 
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READLN (ni) ; 
np := ni*(ni-1) DIV 2; 
T_TABLE_GENERATION; 


{* WRITELN(DATA_FILE, ni:5, np:5, nt-1:5); *) 
NODE_GENERATION; 
PAIR_GENERATION; 
FOR t := 1 TO (nt-1) DO 
WRITELN (DATA_FILE, T_TABLE[t].CAPACITY:6, 
T_TABLE[t])-.COST: 5); 


CLOSE (DATA_FILE) ; 


END; 











C. SOLVING LAGRANGIAN PROBLEM PROCEDURE 


(RRR REE EREREREEEEKREEEREKEREEEEEEEEEREEEERKEEREEEERKERERE ) 


{* TITLE : Generate source-destination pair *} 
{* AUTHOR : Hwang, Yong Goo *) 
{* DATE : 18 July 1990 *) 
{* REVISE : 18 July 1990 *} 


{* DESCRIPTION : This function provides a node pair *) 
{ RARER REKEKEKEKEEREREEEKEREEEEREREREKEEREKEEREREEKERERE | 


FUNCTION p_(i,j: NODE): PAIR; 
VAR 


NODE; 


a: 
: INTEGER; 


oO, 
Pp 
BEGIN 


IF i < j THEN 


2=ptl; 
UNTIL (ORGN[(p] = 0) and (DSTN(p) = d); 
P_:=p; 


END; 


{ RREREKKREEEKEREKEEREEEKEEEREREEEREKEKEREEEEREKEKEKEKKES | 


{* TITLE : Get data and Make route table *) 
{* AUTHOR : Hwang, Yong Goo *) 
{* DATE : 18 July 1990 *) 
{* REVISE : 12 Oct. 1990 *) 








{* DESCRIPTION : This procedure make route table *) 


{* from ‘net.dat’ file *} 
{ RHR EKRKKEKKEKEKKEKKEEEEEEEEKEKEKRERKEREREKERKEEKKERER } 


PROCEDURE GET_DATA; 


VAR 
p: PAIR; 
t: CNFG; 
datf: text; 


(REREKEREREKEKEEKEREREREKE } 
PROCEDURE BUILD_ROUTE_TABLE; 


TYPE 
rte = RECORD 
11,12,13,14: PAIR; 
das: REAL; 
END; 


ptr = “rte; 


VAR 
i,j,k: NODE; 
p,1k1,1k2,1k3,1k4: PAIR; 
c: INTEGER; 
dists: REAL; 
ptrs: ARRAY({1..max_nc] of ptr; 
aptr: ptr; 
nset: SET OF NODE; 
LONGER: BOOLEAN; 


{ RREREKERKERKREKEREKEREKEREEERERE ) 
PROCEDURE REVISE_PTRS (k: INTEGER) ; 


VAR 
c: INTEGER; 


BEGIN 


IF k < ne THEN 
IF dists < ptrs[{k]*.ds THEN 
BEGIN 
dispose(ptrs[{nc-1]); 
FOR c:= nc-1 downto k+1 do 
ptrs[c):=ptrs[c-1]; 
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new(aptr) ; 
with aptr’ do 
BEGIN 


11l:=1k1; 12:=1k2; 13:=1k3; 14:=1k4; 


das:=dists; 
END; 
ptrs[(k]:=aptr; 
END 
ELSE 
BEGIN 
REVISE_PTRS (k+1); 
END 
ELSE 
LONGER: =true; 


END; 


BEGIN {* main of build route table 
{* create number of nc shortest routes 
{* - with the shortest being the direct link. 
WRITELN(’ Wait a moment !! ’); 
FOR p:=1 to np do 
BEGIN 
FOR c:=1 to nc-1 do 
BEGIN 
new(aptr) ; 
aptr’.ds:=VLARGE*10; 
ptrs(c}:=aptr; 
END; 
nset:=[1..niJ; 
nset:=nset-(ORGN(p]]; 
nset:=nset-(DSTN[p]]; 
FOR i:=1 to ni do 
IF i in nset THEN 
BEGIN 
1k1:=p_(ORGN[p],i); 
1k2:=p_(i,DSTN[p]) ; 
1k3:=0; 
1k4:=0; 
dists:=DISTANCE[1k1]+DISTANCE[1k2]; 
LONGER: =false; 
REVISE_PTRS (1); 
IF LONGER = false THEN 
FOR j:=1 to ni do 


44 





IF j in nset-{i] THEN 
BEGIN 
1k2:=p_(i,j)? 
1k3:=p_(j,DSTN(p]}) ; 
1k4:=0; 
dists:=DISTANCE([1k1]+DISTANCE[1k2] 
+DISTANCE[1k3]; 
REVISE_PTRS (1); 
IF LONGER = false THEN 
FOR k:=1 to ni do 
IF k in (nset-[i])~-[j] THEN 
BEGIN 
1k3:=p_(j,k); 
1k4:=p_(k,DSTN[p]}) ; 
dists:=DISTANCE[1k1)+DISTANCE([1k2] 
+DISTANCE(1k3 ]+DISTANCE(1k4]; 
REVISE_PTRS (1); 
END; 
END; 
END; 
ROUTE[ (p-1) *nc,1}:=p 
ROUTE[ (p-1) *nc, 2}: =6 
ROUTE[ (p-1) *nc,3)}:=0; 
ROUTE[ (p-1) *nc, 4] :=0 
FOR c:=1 to nc-1 do 
BEGIN 
ROUTE[ (p-1) *nc+c,1)]:=ptrs[(c])*.11; 
ROUTE[ (p-1) *nc+c,2):=ptrs[c})*.12; 
ROUTE[ (p-1) *nc+c,3}:=ptrs[(c}*.13;7 
ROUTE[ (p-1) *nce+c,4]:=ptrs[(c]*.14; 
END; 


END; (* FOR p *} 


FOR c:=1 to nc-1 do 
BEGIN 
dispose (ptrs[c]); 
ptrs(cJ:=nil; 
END; 


END; {* end of build route table *} 


BEGIN {* main of get data *} 





ASSIGN (datf, F_NAME) ; 
RESET (datf) ; 


FOR p:=1 to np do 
READLN (datf,ORGN[p] ,DSTN[p],DISTANCE[p],TRAFFIC[p}); 


WRITELN(’ Enter the No. of Candidate routes for each Pair :’); 
READIN (nc) ; 
CAPA[0]:=0; 
RATE[0}:=0; 


FOR t:=1 to nt do 
READLN (datf,CAPA[t],RATE[t]); 


BUILD ROUTE_TABLE; 


END; 


{ RHRRERKKKEK KEKE KEK KRERE KEKE EEK EKKKEEKREKKEEKKKKEREKKEKE | 


{* TITLE : Find initial feasible solution *} 
{* AUTHOR : Hwang, Yong Goo *} 
{* DATE : 2 Aug. 1990 *} 
{* REVISE  : 15 Aug. 1990 *} 
{* DESCRIPTION : This procedure gives us initial *} 
{* feasible solution for problem. This *)} 
{* solution may be fully connected network. *) 


{ RREKKKKEKKKEKKEEKEREKREREEKKEEKEREEREEKEEKEKEEKEKEKEKEKKEE } 


PROCEDURE FIND_INITIAL_FEASIBLE SOLUTION; 
VAR 

Pp: PAIR; 

t: CNFG; 
BEGIN 

pobj :=0; 


FOR p:=1 to np do 
BEGIN 
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:=0; 
IF TRAFFIC({p] > 0 THEN 
BEGIN 
REPEAT 
2=t+1; 
UNTIL TRAFFIC(p] <= CAPA[t]; 
END; 
px[p]:=(p-1) *nc; 
pZ2[p):=t; 
pobj :=pobj+DISTANCE[p])*RATE[(t]; 
END; 
UPPBND:=pobj ; 


END; 


{ REHRKRHK KEKE KEKE KEKE KEKE KEKKEKEK KEKE REE REKKEKEREKEREE | 


(* TITLE : Initiate subgradient parameters *} 
{* AUTHOR : Hwang, Yong Goo *)} 
{* DATE : 13 June 1990 *} 
{* REVISE : 15 Nov. 1990 *) 
{* DESCRIPTION : At first processing, this *} 
{* procedure initiates sub. parameters *} 
{* autometically, and then when user reenter *} 
{* the system, he can modify them. *} 


{ Keka aaa KKH KEE KKK KKK KKK KKKKEKEKKEEKKEKKKEKEKKEKE } 


PROCEDURE PREP_FOR_SUBGRADIENT_ITERATION; 


VAR 
p : PAIR; 
ANS : CHAR; 


BEGIN 


IF FRAG] THEN 
BEGIN 

LAMBDA:= 0.5; 
IMPLIM:= 25; 
ITRLIM:=100; 
ERRORLIM: =0.0050; 
LAMBDALIM:=0.0001; 
ITRCTR: =0; 
IMPCTR:=0; 
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LOWBND: =0; 


BEGIN 
WRITELN(’Do you want to change some subgradient parameters? ‘); 
READLN (ANS) ; 
IF (ANS = ‘’y’) OR (ANS = ‘’Y’) THEN 
BEGIN 
WRITELN(’ LAMBDA : ’, LAMBDA :4:2); 
WRITELN(’ Enter the LAMBDA? ’); 
READLN (LAMBDA) ; 
WRITELN(’ IMPROVEMENT COUNTER LIMIT : ’, IMPLIM :4); 
WRITELN(’ Enter the IMPROVEMENT COUNTER LIMIT? ’); 
READLN ( IMPLIM) ; 
WRITELN(’ LAMBDA LIMIT : ’, LAMBDALIM :4:6); 
WRITELN(’ Enter the LAMBDA LIMIT? ‘); 
READLN (LAMBDALIM) ; 


WRITELN(’ ITERATION COUNTER LIMIT : ’, ITRLIM :4); 
WRITELN(’ Enter the ITERATION COUNTER LIMIT? ’); 
READLN (ITRLIM) ; 
END 
ELSE 
BEGIN 


LAMBDA:= 0.5; 
IMPLIM:= 25; 
ITRLIM:=100; 
LAMBDALIM:=0.0001; 
END; 
ERRORLIM:=0.0050; 
ITRCTR:=0; 
IMPCTR:=0; 
LOWBND:=0; 
FOR p:=0 to np do 
ALPHA[p] :=0; 
END; 


END; 


{ RERERKKKKEKEREREEEREKREREEEEKEKEKEEREEEEEEEERERERERERES } 
{* TITLE : Solve the Lagrangian dual *) 
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{* AUTHOR : Hwang, Yong Goo *} 
{* DATE : 21 Sep. 1990 *) 
{* REVISE : 11 Nov. 1990 *} 
{* DESCRIPTION : This procedure computes *) 
{* a Lagrangian dual. *)} 


{ ERREKRKKKKKEKKKKEKKKEEEKEKEKEKEREEKKEEEKEEKREEEKEEEEKKEEEKEEE ) 


PROCEDURE SOLVE_LAGRANGIAN_ DUAL; 


VAR 
Pp: PAIR; 
t: CNFG; 
r: INTEGER; 
sum,min: REAL; 


BEGIN 


dobj :=0; 
FOR p:=1 to np do 
BEGIN 
min: =VLARGE; 
FOR r:=(p-1)*nc to p*nc-1 do 
BEGIN 
sum:=ALPHA[ROUTE[r,1]]+ALPHA[ROUTE[r,2)}] 
+ALPHA[ROUTE[r,3]]+ALPHA[ROUTE(r,4]]; 
IF sum < min THEN 
BEGIN 
min:=sum; 
ax[(p]:=r; 
END; 
END; 
dobj :=dobj+TRAFFIC[p)*min; 
END; 
FOR p:=1 to np do 
BEGIN 
min:=0; 
az(pj]:=0; 
FOR t:=1 to nt do 
BEGIN 
IF DISTANCE[p]*RATE[t]-ALPHA[p]*CAPA[t] < min THEN 
BEGIN 
min: =DISTANCE[p]*RATE[t]-ALPHA[p) *CAPA[t]; 
az(pj:=t; 
END; 
END; 
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aobj :=dobj+min; 
END; 
END; 


{ ERHERAKEKEKEKEKEKEKEEKEKEEREKEEKERERREREREKEKEEKREREEEKES ) 


{* TITLE : Make feasible solution *) 
{* AUTHOR : Hwang, Yong Goo *)} 
(* DATE : 15 Sep. 1990 *} 
{* REVISE : 15 Nov. 1990 *} 
{* DESCRIPTION : This procedure will give us *} 
{* a feasible solution using infeasible sol- x} 
{* ution. If this solution is lower than *) 
{* upper bound, then it replaces upper bound. *} 


{ RHRKKKKKEKKEKEKKEKEREKKEKKEKEEKEKKKKKEKKEKERKEKEREEKEKEKE ) 


PROCEDURE PERTURB_FOR_FEASIBLE_SOLUTIONS; 


VAR 
p: PAIR; 
t: CNFG; 
obj: REAL; 
flag: BOOLEAN; 


BEGIN 


obj :=0; 
FOR p:=1 to np do 
BEGIN 
flag:=false; 
t:=0; 
IF SUBGRAD[p)+CAPA[daz[p]] > 0 THEN 
BEGIN 
REPEAT 
IF t+1l <= nt THEN 
ts=t+1 
ELSE 
flag:=true; 
UNTIL flag or (SUBGRAD[p]+CAPA[dz[{p]] <= CAPA[t]); 
END; 
IF flag THEN 











obj :=obj}+DISTANCE[p)*RATE[t]; 
END; 
END; 
IF obj < UPPBND THEN 
BEGIN 
UPPBND: =obj ; 
pz:=daz; 
px: =ax; 
END; 


END; 


{ REKKKKKKKKKKKKREKKKKEKEKKKKEKKEKEKEKEEEKEKEKEKEKKEKEKEKEE ) 


{* TITLE : Test subgradient parameters *)} 
{* AUTHOR : Hwang, Yong Goo *} 
{* DATE : 15 Sep. 1990 *) 
{* REVISE : 18 Nov. 1990 *) 
{* DESCRIPTION : This procedure checks that *)} 
{* terminated condition is satisfied. *} 


{ RRR AKEKKKKKKKEKEKKKEKEKKKKKKEEKEEEKEKKEKEKEEKEKKEKKEKKER | 


PROCEDURE CHECK_IMPRO_COUNTER; 
BEGIN 


IF IMPCTR > IMPLIM THEN 
BEGIN 
LAMBDA: =LAMBDA/2.0; 
IMPCTR:=1; 
END; 


END; 


( RRRRERREREKEEKEREEREEEKKE } 
PROCEDURE TEST_TERMINATION; 


BEGIN 


IF (LOWBND > 0) and (UPPBND-LOWBND > 0) 
and ((UPPBND-LOWBND) /LOWBND < ERRORLIM) THEN 
TERMINATED: =true 
ELSE 


IF LAMBDA < LAMBDALIM THEN 











TERMINATED: =true 
ELSE 
IF ITRCTR >= ITRLIM THEN 
TERMINATED: =true 
ELSE 
TERMINATED: =false; 


END; 


{ REHRERAEKEKEKEEKEKEKEKEREREREREREEEEEEEEEREREEEEREKEKEEKE | 


{* TITLE : Update Lagrange multipliers *)} 
{* AUTHOR : Hwang, Yong Goo *} 
{* DATE : 12 Aug. 1990 *) 
{* REVISE : 13 Nov. 1990 *)} 
{* DESCRIPTION : This procedure change Lagrange *} 
{* multipliers. *) 


{ RRR KEEEKEEEKEKEKEEKEKEREKEKEEEEEEEKEEEEEEKEREKKEKESE } 


PROCEDURE UPDATE_LAGRANGE_MULTIPLIERS; 


VAR 
p: PAIR; 


BEGIN 


STEPSIZE:=LAMBDA* (UPPBND-dobj ) /NORM; 
FOR p:=1 to np do 
BEGIN 
ALPHA([p) :=ALPHA[p]+STEPSIZE*SUBGRAD([p}; 
IF ALPHA[p] < le~4 THEN 
ALPHA[p] :=0; 
IF ALPHA[p] >= VLARGE THEN 
BEGIN 
WRITELN(’WARNING: MULTIPLIER TOO LARGE’) ; 
WRITELN (dobj :10:2,NORM:10:0,STEPSIZE:10:2,SUBGRAD[p]:10:2); 
READLN ; 
END; 
END; 


END; 
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{ RRRAKKKKKRKKKKKKKKEEKEKKKKKKEKKKKKKKKKKKKKKKKKKKKKKKKKKS ) 


TITLE : Solve problem 
AUTHOR : Hwang, Yong Goo 
DATE : 7 Oct. 1990 
REVISE : 16 Nov. 1990 


DESCRIPTION : This procedure make a best 
feasible solution using above all subproc- 
edure. 


{ RRKKKKKKRRAKEKERKEKKKEKEKKEEKEEKEREKEKEEKKKKKEKKKKKEKKKKKE } 


PROCEDURE SOLVE_PROBLEM; 


BEGIN 


REPEAT ; 


ITRCTR: =ITRCTR+1; 
IMPCTR: =IMPCTR+1; 


CHECK_IMPRO_COUNTER; 


SOLVE_LAGRANGIAN_ DUAL; 


WRITELN(ITRCTR:4,’ ‘’,dobj:10:2,’ ‘’,UPPBND:10:2,’ 


IF dob} > LOWBND THEN 
BEGIN 
IMPCTR:=0; 
LOWBND: =dobj; 
BEST_ALPHA:=ALPHA; 
END; 


SUBGRAD[0]:=0; 
FOR p:=1 to np do 
SUBGRAD[p] :=CAPA[dz[p])]*(-1); 
FOR p:=1 to np do 
BEGIN 
SUBGRAD[ROUTE[ax[p],1] 
SUBGRAD[ROUTE[dax[p]},2] 
SUBGRAD[ROUTE[(dx[p],3] 
SUBGRAD[ROUTE[dx{p] ,4] 
END; 
FOR p:=1 to np do 
IF abs(SUBGRAD[(p}) < le-4 THEN 
SUBGRAD[p]:=0; 


—_ 
= 
or 
¢ =, 
c= 
e 


ed heed Od Od 
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+ LAMBDA: 10:4); 


SUBGRAD[ROUTE[dx[p],1]]+TRAFFIC[p]; 
SUBGRAD[ROUTE[dx[{p],2]]+TRAFFIC(p}; 
SUBGRAD[ROUTE[dx[(p],3]]+TRAFFIC[p]; 
SUBGRAD[ROUTE[dx[p],4]]+TRAFFIC[p]; 


NORM: =0; 
FOR p:=1 to np do 

NORM: =NORM+SUBGRAD[p] *SUBGRAD[p] ; 
OPTIMAL: =true; 
FOR p:=1 to np do 

IF SUBGRAD[p] > O THEN 

OPTIMAL: =false; 

IF OPTIMAL THEN 


BEGIN 
UPPBND: =dob3 ; 


PERTURB_FOR_FEASIBLE_SOLUTIONS; 
TEST_TERMINATION; 


IF not TERMINATED THEN 
UPDATE_LAGRANGE MULTIPLIERS; 


UNTIL TERMINATED; 


END; 
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D. REPORT PROCEDURE 


{ ERRAKKEKKKEKEKKRKKKKEREKKEKEKEREEREEKRERREKERKEKREKEKEEKKEKE } 


{* TITLE : Make report *) 
{* AUTHOR : Hwang, Yong Goo *} 
{* DATE : 3 Oct. 1990 *) 
{* REVISE : 21 Nov. 1990 *)} 


{* DESCRIPTION : This procedure provides a report. *)}) 
{ RHRRAKEKKKKAKKKKEKKKKEKKEEKREKKERREKEKKKKKKEKEKKKKKKKKKKE ) 


PROCEDURE DRAW_NETWORK; 


CONST 
MAXLENGTH = 250; 
MAXHEIGHT = 150; 


VAR 
CENTER_X, 
CENTER_Y : INTEGER; 
i, j : INTEGER; 
TEMP_FT_1, TEMP_T_1, TEMP_T_3 : INTEGER; 
inkey : INTEGER; 
ckey : char; 
DRIVERVAR, MODEVAR : INTEGER; 


(RERRERKERREREKRREREREREEREREREEEKEREKE ) 
PROCEDURE GRAPHTITLE (INTITLE : STRING) ; 


BEGIN 
SETTEXTJUSTIFY (CENTERTEXT, TOPTEXT) ; 
SETTEXTSTYLE (TRIPLEXFONT, HORIZDIR, 4); 
OUTTEXTXY (GETMAXX DIV 2, 1, INTITLE); 
END; {* END of graphtitle *)} 
{ KERKEKKKEKEKKKKKKKKKKE } 


PROCEDURE GRAPHCONTINUE; 


CONST 
messagel = ‘CONTINUE? (y/n)/’; 
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MESSAGE2 = ‘Enter the PAIR number!’; 
BEGIN 


SETTEXTJUSTIFY (CENTERTEXT, TOPTEXT) ; 
SETTEXTSTYLE (GOTHICFONT, HORIZDIR, 3); 


OUTTEXTXY (GETMAXX DIV 2, GETMAXY - 40, MESSAGE1); 


READLN (ckey) 7 
IF (ckey <> ’N’) AND (ckey <> ‘n’) THEN 
BEGIN 
SETTEXTJUSTIFY (CENTERTEXT, TOPTEXT) ; 
SETTEXTSTYLE (GOTHICFONT, HORIZDIR, 3); 


OUTTEXTXY (GETMAXX DIV 2, GETMAXY - 50, MESSAGE2) ; 


READLN (inkey) ; 
END; 
CLEARDEVICE; 


END; {* END of graphcontinue *) 


( BEKKKEKEEEKEKEKEKERE } 
PROCEDURE DRAW_ROUTE; 


VAR 
i,j 2 INTEGER; 


BEGIN 


LINE (10,10,10,450); 

LINE (10,450,getmaxx, 450) ; 

LINE (center_x,450,center_x,10); 
LINE (getmaxx,10,10,10) ; 

LINE (getmaxx,getmaxy,getmaxx,10) ; 
FOR i := 1 TO ni DO 


CIRCLE((POINT_X[i] DIV 3),(POINT_Y[{i] DIV 3), 2); 


FOR j := 1 TO np DO 
IF (pz{j] <> 0) THEN 
BEGIN 
{* TEMP FT 1 := 
TEMP_T 1 := T_TABLE[{j].T_1; 
TEMP T 3 := T _TABLE[j].T_3; *} 
LINE(POINT_X[(S_D_PAIR{j}}.S] DIV 3, 
POINT Y{S_D_PAIR[J].S] DIV 3, 
POINT _X[S_D_PAIR{j].D] DIV 3, 
POINT_Y{S_D_PAIR[J].D] DIV 3); 


T_TABLE(j].FT_1; 














END; 


LINE (POINT_X{s_d_PAIR[ROUTE[px[inkey],1]].s] div 
POINT _Y({s_d_ PAIR(ROUTE[px[inkey],1]].s] div 
POINT_X[s_a | PAIR[ROUTE[px[inkey], 1})}.d}) div 
POINT_Y[{s_d | PAIR[ROUTE[px[inkey], 1}}.daj} div 


IF ROUTE(px{inkey],2] > 0 THEN 


LINE(POINT_X[{s_d_PAIR[ROUTE(px[inkey],2]].s] 
POINT _Y[s_d_PAIR[ROUTE[(px[inkey],2]}].s] 
POINT_X[s_d_PAIR[ROUTE[px[inkey},2]]}.d] 
POINT _Y[s_d_PAIR[ROUTE{px[inkey],2]}].d] 


IF ROUTE(px{inkey],3] > 0 THEN 


LINE (POINT _X[s_d_PAIR[ROUTE[px[inkey],3]]-.s] 
POINT Y[s_d_PAIR[ROUTE[px[inkey],3]}-s] 
POINT _X[s_d_PAIR[ROUTE[px[inkey},3]].a} 
POINT Y[s_d_PAIR{ROUTE[px[inkey],3]}-d] 


IF ROUTE[px{inkey],4] > 0 THEN 


LINE (POINT_X[s_d_PAIR[ROUTE[px{inkey],4]]-s] 
POINT _Y[{s_d_PAIR(ROUTE(px([inkey],4]].s] 
POINT _X[S_d_PAIR[ROUTE[px[inkey],4]]-d] 
POINT Y[{s_d_PAIR[{ROUTE[px[inkey],4]]-.d] 


END; 


BEGIN {* MAIN PART *} 


inkey :=0; 

DRIVERVAR := detect; 
INITGRAPH(DRIVERVAR, MODEVAR, ‘'’); 
CENTER_X := GETMAXX DIV 2; 
CENTER_Y := GETMAXY DIV 2; 
GRAPHTITLE(’NETWORK DESIGN’) ; 
LINE(10,10,10,450); 

LINE (10,450,550,450); 

LINE (550,450,550,10); 

LINE (550,10,10,10); 

FOR i := 1 TO ni DO 





4 + center_x, 


4, 
4 + center_x, 
4); 
div 4 + center _x, 
div 
div center_x, 
div 
div center_x, 
div 
div center x, 
div 
div center _x, 
div 
div center_x, 
div 


CIRCLE ((POINT_X[i]) DIV 2),(POINT_Y[i] DIV 2), 


FOR 3 := 1 TO np DO 


BEGIN 
IF (pz(j] <> 0) THEN 
BEGIN 
{* TEMP_FT_1 := T_TABLE[j}].FT_1; 
TEMP_T 1 := T_TABLE[}].T_1; 
TEMP T 3 := T_TABLE(j].T_3; *} 
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LINE (POINT _X(S_D_PAIR(j}.S) DIV 2, 
POINT Y[(S_D PAIR[J].S] DIV 2, 
POINT X[S_D_PAIR(j].D] DIV 2, 
POINT Y[S_D PAIR[J}.D] DIV 2); 
END; 
END; 
REPEAT 
graphcontinue; 
IF (inkey <> 0) THEN 
DRAW_ROUTE; 
UNTIL (ckey = ’N’) or (ckey = ’n’); 
CLOSEGRAPH ; 


END; 


{ RRERK AKER REKEKEREKKERERKEKEEKEREREE } 
PROCEDURE REPORT_BEST_FEASIBLE_ SOLUTION; 


VAR 
p: PAIR; 
af: text; 
i: INTEGER; 


BEGIN 


ASSIGN (df, FILE_NAME) ; 
REWRITE (af) ; 
WRITELN (df); 
WRITELN (df, ’relative error bound: ’, 
( (UPPBND~-LOWBND) /UPPBND) *100:7:4,/’ %/’); 
WRITELN (af, ’LAMBDA: ’, LAMBDA) ; 
WRITELN (df) ; 
WRITELN(df,’The best feasible solution: ’,UPPBND:10:2); 
FOR p:=1 to np do 
BEGIN 
WRITE(df,p:2,’ [’,ORGN(p]:2,’-’,DSTN[p]:2,’] : ’ 
,ROUTE[px[p],1]:4); 
IF ROUTE(px[(p],2] > 0 THEN 
WRITE (df,ROUTE[px[p],2]:4); 
IF ROUTE(px[p],3] > 0 THEN 
WRITE (af, ROUTE[px[p},3]:4); 
IF ROUTE[px{p],4] > 0 THEN 
WRITE (df, ROUTE[px[p],4]:4); 
WRITELN (df) ; 
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END; 
i:=0; 
FOR p:=1 to np do 
BEGIN 
IF i < 4 THEN 
BEGIN 
WRITE (Qf,’[',p:3,’] ‘,pz{p}:4); 
i:=i+1; 
END 
ELSE 
BEGIN 
WRITE (dGf,’[’,p:3,’] ’,pz[{p]:4); 
WRITELN (af); 
i:=0; 
END; 
END; 
WRITELN (af) ; 
close(af); 
DRAW_NETWORK; 


END; 


{ SRERKEKKKKKKKEKKKKKEKEKKKKKKKKEKEKKKEKE } 
{* main body *) 
{ RHRAEKKKKRAKKKEKKKKKKKKKKKKKKERKEKEKEKE | 
BEGIN 
FRAG1 := TRUE; 
GENERATE_DATA; 
REPEAT 
GET_DATA; 
FIND _INITIAL_ FEASIBLE SOLUTION; 
PREP_FOR_SUBGRADIENT_ITERATION; 
SOLVE_PROBLEM; 


REPORT_BEST_FEASIBLE_SOLUTION; 
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WRITELN(’ Try again? : (y/n) ’)? 
READLN (ANS1) ; 


UNTIL (ANS1 = ’N’) OR (ANS1 = /n’); 


END. {*** main PROCEDURE ***)} 
{BHR EREREKEKEREKEREKEKKEEREREKEKEREREKREREEEEERE } 
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APPENDIX B: SAMPLE INPUT DATA FILE 


A. PAIR INFORMATION 


pairs 
sour. dest. distance traffic 


1 2 395 7 
1 3 339 98 
1 4 156 0 
1 5 261 64 
1 6 724 43 
1 7 516 23 
2 3 556 0 
2 4 448 38 
2 5 592 27 
2 6 674 0 
2 7 502 93 
3 4 491 62 
3 5 172 90 
3 6 469 28 
3 7 295 34 
4 ~ 385 53 
4 6 877 0 
4 7 668 87 
5 6 640 23 
5 7 459 0 
6 7 209 6 
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B. CIRCUIT INFORMATION 


CIRCUIT 
capa. cost 
6 10 
12 15 
18 19 
24 21 
30 31 
36 36 
42 40 
48 42 
54 52 
60 57 
66 61 
72 63 
78 73 
84 78 
90 82 
96 84 
102 94 
108 99 
114 103 
120 105 
126 115 
132 120 
138 124 
144 126 
150 136 
156 141 
162 145 
168 147 
174 157 
180 162 
186 166 
192 168 
198 178 
204 183 
210 187 
216 189 
222 199 
228 204 
234 208 
672 219 
678 229 
684 234 
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APPENDIX C: SAMPLE OUTPUT DATA FILE 


relative error bound: 63.5411 % 
LAMBDA: 2.5000000000E-01 


The best feasible solution: 274700.00 


1 [ i= 2} 3 1 
2 [1-3]: 4 13 
3 [ 1- 4} 3 
4 [1-5] 4 
5 [ 1- 6] 4 13 14 
6 [ 1- 7] 4 13 15 
7 ( 2- 3] ° 
8 [ 2- 4]: 1 3 
9 [ 2- 5): 1 4 
10 [ 2- 6} : 10 
11 [ 2-7): 112 
12 [ 3-4) : 13 4 3 
13 [ 3-5) : 13 
144 [3-6]: 14 
15 [3-7] : 15 
16 [ 4- 5): 3. 4 
17 [ 4- 6} : 17 
18 [ 4-7): 3. 4 13 15 
19 [ 5- 6] : 13 14 
20 [ 5-7) : 13 15 
21 ( 6-7]: 212 
aes ees 1 Ge Of 3] 40[ 4} 40[ 5] 0 
[56] Of 7] Of 8] of 9] O[ 10) 0 
{11} 16[ 12} Of 13} 40[ 14) 16[ 15] 24 
[ 16] Of 17] O[ 18] O[{ 19} O[ 20] 0 
{ 22) 1 





1]. 
12. 


LIST OF REFERENCES 


P. Afentakis and B. Gavish, ”Optimal lot-sizing algorithms for complex product 
structures”, Grad. School Management, University Rochester, Rochester, NY, 
1983. 


. D.P. Bertsekas, ”A class of optimal routing algorithms for communications net- 


works”, in Proc. 1980 Int. Conf. Circuits Comput., Atlanta, GA, Nov. 1980 


. MLL. Fisher, ”Lagrangean relaxation method for solving integer programming 


problems”, Management Science, vol. 27, pp. 1-18, Jan. 1981. 


. M. Held, P. Wolfe, and H. Crowder, ”Variation of subgradient optimization”, 


Math. Program., vol. 6, pp. 62-88, 1974. 


. B. Gavish and S.L. Hantler, ”An algorithm for optimal route selection in SNA 


networks”, IEEE Trans. on Comm., vol.com-31, pp. 1154-1160, Oct. 1983. 


. M. Karnaugh, ”A New Class of Algorithms for Multipoint Network Optimiza- 


tion”, IEEE Transactions on Communications, vol. COM-24, pp. 500-505, May 
1976. 


. L.G. Mitten, ”Branch-and bound methods: General formulation and proper- 


ties”, Operations Research, vol. 18, pp 24-34, 1970. 


. E. Rosenberg, ”A Nonlinear Programming Heuristic for Computing Optimal 


Link Capacities in a Multi-hour Alternate Routing Communications network”, 
Operations Research, vei. 35, No. 3, pp. 354-364, May-June 1987. 


. D.N. Lee, K.T. Medhi and J.L. Strand, "Solving Large Telecommunications 


Network Loading Problems”, AT T Technical Journal, pp. 48-56, May-June 
1989. 


. W.J. Barksdale, Practical Computer Data Communications, Plenum Press, New 


York, 1986. 
W.A. Flanagan, The Guide to T-1 Networking, Telecom Library Inc., 1988. 


R. Sharma, ”T-1 Network Design and Planning Made Easier”, Data Commu- 
nications, pp. 199-207, Sept. 1988. 


66 


13. 


14. 


' 15. 


16. 
17. 





S. Fleming, "Get Ready for T-3 Networking”, Data Communications, pp. 82-94, 
Sept. 1989. 


Udi Manber, Introduction to Algorithms, Addison-Wesley Publishing Co., July 
1988. 


David R. Anderson, Dennis J. Sweeney and Thomas A. Williams, An Introduc- 
tion to Management Science, West Publishing Co., 1988. 


Andrew S. Tanenbaum, Computer Networks, Prentice-Hall Inc., 1988. 
Uyless Black, Data Networks, Prentice-Hall Inc., 1989. 


67 





i) 


INITIAL DISTRIBUTION LIST 


. Defense Technical Information Center 


Cameron Station 
Alexandria, Virginia 22304-6145 


. Library, Code 52 


Naval Postgraduate School 
Monterey, California 93943 


. Department Chairman, Code CS 


Department of Computer Science 
Naval Postgraduate School 
Monterey, California 93943 


. Prof. Myung W. Suh, Code ASSU 


Department of Administrative Sciences 
Naval Postgraduate School 
Monterey, California 93943 


. Prof. Richard W. Hamming, Code CSHG 


Department of Computer Science 
Naval Postgraduate School 
Monterey, California 93943 


. Prof. Mantak Shing, Code CSSH 


Department of Computer Science 
Naval Postgraduate School 
Monterey, California 93943 


. Prof. Dong S. Kim, Code MEKM 


Department of Mechanical Engineering 
Naval Postgraduate School 
Monterey, California 93943 


. Prof. Kyung C. Kim, Code 


Department of Computer Science 
Naval Postgraduate School 
Monterey, California 93943 


68 








10. 


11. 


12. 


13. 


14. 


. Army Central Library 
Army Headquarter, Bunam-Ri, 

Duma-Myun, Nonsan-Gun, Chungcheong-Nam-Do, 
Republic of Korea, 320-919 


Hwang; Yong Goo 


734-3, HwangGeum-Dong, SuSeong-Gu, Teagu, 
Seoul, Republic of Korea, 706-040 


Song, 1 Youn 
SMC 2302 NPGS 
Monterey, California 93943 


Chung, Jae Du 
SMC 1504 NPGS 
Monterey, California 93943 


Park, Hyun Kyoo 
SMC 2820 NPGS 
Monterey, California 93943 


Choi, Nag Jung 
SMC 2853 NPGS 
Monterey, California 93943 








69 





